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ABSTRACT 

The inverse Compton scattering of the cosmic microwave background (CMB) 
radiation with electrons in the intracluster medium which has a temperature 
gradient, was examined by the third-order perturbation theory of the Compton 
scattering. A new type of the spectrum distortion of the CMB was found and 
named as gradient T Sunyaev-Zel'dovich effect (gradT SZE). The spectrum has 
an universal shape. The spectrum crosses over zero at 326GHz. The sign of the 
spectrum depends on the relative direction of the line-of-sight to the direction 
of the temperature gradient. This unique spectrum shape can be used to detect 
the gradT SZE signal by broad-band or multi-frequency observations of the SZE. 
The amplitude of the spectrum distortion does not depend on the electron density 
and is proportional to the heat conductivity. Therefore, the gradT SZE provides 
an unique opportunity to measure thermally nonequilibrium electron momentum 
distribution function when the ICM has a temperature gradient and the heat 
conductivity in the ICM. However, the expected amplitude of the signal is very 
small. The modifications to the thermal SZE spectrum due to variety of known 
effects, such as relativistic correction etc., can become problematic when using 
multi-frequency separation techniques to detect the gradT SZE signal. 

Subject headings: galaxies: clusters: general — magnetic fields — conduction— 
Compton scattering — plasma-mm and sub-mm observations 
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1. Introduction 

The recent dramatic progresses of the X-ray observations have been unveiling that the 
temperature distribution of the intracluster medium (hereafter ICM) is far from isother- 
mal. It has turned out that the heat conductivity is one of the key parameters controlling 
the evolution of the ICM(Ikebe et al. 1999; Markevitch et al. 2000; Fabian et al. 2001; 
Vikhlinin et al. 2001; Markevitch et al. 2003; Zakamska & Narayan 2003). The global 
temperature gradient have been found in some cooling core clusters(Kaastra et al. 2004). 
Zakamska & Narayan (2003) performed hydrostatic model fitting of the observed tempera- 
ture and density profile of the ICM. They assumed the energy balance between the radiative 
cooling and the conductive heating. They showed that the heat conductivities should be one- 
third of the Spitzer value in majority of the clusters. Shock waves(Markevitch et al. 2002; 
Fujita et al. 2004) and cold fronts(Markevitch et al. 2000; Vikhlinin et al. 2001) found in 
merging clusters are showing the temperature jumps across the discontinuities. The hot 
bubbles caused by the interaction with the activity of the central galaxies have been found 
in many clusters(McNamara et al. 2000; Fabian et al. 2000; Blanton et al. 2001). The heat 
conductivities across the cold fronts and the surface of the hot bubbles have to be reduced in 
many order of magnitudes to maintain the structures(Ettori & Fabian 2000). The large scale 
fluctuation of the temperature distributions have been found in several clusters(Watanabe et 
al, 1999; Shibata et al. 2001; Markevitch et al. 2003). To maintain the observed temperature 
fluctuations at least for a dynamical time scale of clusters, the order of magnitude reduction 
of the heat conductivities are required(Markevitch et al. 2003). All the above mentioned 
facts show that a direct method for measuring the heat conductivities may provide crucial 
informations to study the evolution of the ICM. 

The recent theoretical progresses have been unveiling that a microscopic instability of the 
plasma which has a temperature gradient (Ramani & Laval 1978), could play important roles 
on the suppression of the heat conductivities(Levinson & Eichler 1992; Hattori & Umetsu 2000; 
Okabe & Hattori 2003) and the origin of cluster magnetic fields(Okabe & Hattori 2003; Ok- 
abe & Hattori 2004). The heat conduction is the electron thermal energy flow from hot 
to cold regions. Since it is the third order moment of the electron momentum distribu- 
tion function, the electron momentum distribution function must have deviation from the 
Maxwell-Boltzmann distribution which is the even function of the momentum. Therefore, the 
electron momentum distribution function of the plasma with a temperature gradient is de- 
scribed by the sum of the Maxwell-Boltzmann distribution, g m (q) = (27r) 3 n e (27rm e /csT e )~ 3 / 2 
cxp 2m 9 k B T ) ' an d the non Maxwellian part, Ag(q), where q is a momentum of electron, 
and n e , T e and m e are electron number density, temperature and mass, respectively. The 
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heat current density, nVT e , is described by the Ag(q) as 

where f K is the heat conductivity in the Spitzer value unit and k$ p is the Spitzer heat conduc- 

/, r \l/2 

tivity. The Spitzer heat conductivity is given by ks p = 1.31X e n e kB ( -^- s - ) where A e is the 
electron Coulomb mean free path(Spitzer 1956; Sarazin 1988). The deviation of the electron 
momentum distribution from the thermal equilibrium distribution(Okabe & Hattori 2003) 
in the plasma with a temperature gradient is described by 

* / \ , m q- VT e / 4g 2 \ 

where q t h = V^ m ekBT e is the electron thermal momentum. The plasma instability is derived 
by this local deviation of the electron momentum distribution from the thermal equilibrium 
distribution. It has been shown that the instability may provide an origin of the magnetic 
fields in cluster of galaxies (Okabe & Hattori 2003). Therefore, the measurement of the func- 
tional form of the non-equilibrium momentum distribution function is crucial for confirming 
the above mentioned theoretical models on the suppression of the heat conductivities and 
generation of the magnetic fields. 

In this paper, a spectrum distortion of the cosmic microwave background (hereafter 
CMB) due to the inverse Compton scattering with electrons which support heat conduction, 
is examined by solving the Boltzmann equation for Compton scattering to third-order in the 
electron momentum. The result provides a direct method for measuring the heat conductiv- 
ities and the functional form of the electron momentum distribution functions in the plasma 
with a temperature gradient. 



2. The Boltzmann equation for third-order Compton scattering 

The spectrum distortion of the cosmic microwave background (CMB) due to the Comp- 
ton scattering from electrons with momentum distribution function of g(q) is given by the 
Boltzmann equation(Hu et al. 1994) as 

^5T~ = ^J D * D 4Dp\2n)WHp + q-$-q>)\M\* 

x {<7(x, q')/(x, p') [1 + /(x, p)] - 0(x, q)/(x, P ) [1 + /(x, p')]}, (2) 

where /(x, p) is the photon occupation number with a three- momentum p, p and p' are the 
three-momentum for photons of four-momentum p and p', q and q' are the three-momentum 
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for electrons of four- momentum q and q', E(p) is the energy of the photon with a momentum 
p, \M\ 2 is the Lorentz-invariant matrix element for Compton scattering(Hu et al. 1994), and 

n d3q 

q W2E(q) 

is the Lorentz-invariant phase space volume element of electrons with an energy of E(q) and 
a momentum q. The left hand side of the equation is the total derivative of the photon 
distribution function in the phase space written as, 

Df df_ d£dx^ dj^dn^ Of dEjp) 
~Dt ~ ~dt + dx i ~dt + dn~dt Jr dE(p) dt ' 

where Hi is the i-th component of the direction cosines, n, for a photon of four-momentum 
p. The first two terms in the right hand side is a derivative along the photon path. The 
third term represents an effect of gravitational deflect of the photon path due to the clus- 
ter gravitational lensing effect. Since the deflection angle due to the cluster gravitational 
lensing effect is at most 30", the change of the path length due to the lensing effect is at 
most O.Olpc and negligibly small compared with the cluster size of ~lMpc. The last term 
represents the cosmological redshift and the energy change of the photons due to the grav- 
itational scattering by the tangential peculiar motion of the cluster (Birkinshaw et al. 1984; 
Aso et al. 2002). The gravitational scattering causes the CMB temperature fluctuation of 
ST/T ~ 10~ 7 (Aso et al. 2002). The effect is order of magnitude smaller than the Compton 
scattering effect concerned in this paper as shown in the next section. Therefore, the terms 
representing the gravitational lensing effect and the gravitational scattering effect are ne- 
glected in the course of the following calculations. Further, the cosmological redshift effect 
is neglected for simplicity since the effect is trivial(Hu et al. 1994). 

The right hand side of the Boltzmann equation is the collision term. The first and 
the second terms are the contributions from scattering into and out of the momentum state 
p including the stimulated emission effects, respectively. In the rest frame of the incident 
electron, the matrix element for Compton scattering summed over polarization is given by, 

|M| 2 = 2(4tt)V 

where e is the absolute value of the charge of the electron. The variables with ~ denotes 
those in the electron rest frame, otherwise those in the laboratory frame. 

As an example, the case when the electron velocity distribution function is the Maxwell- 
Boltzmann and the thermal energy contained in electron system is much less than the electron 
rest mass energy, is considered. By expanding the collision term of the Boltzmann equation 
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for isotropically distributed incident photons up to the second-order in electron momentum, 
the Kompaneets equation is deduced(Hu et al. 1994; Kompaneets 1957). The spectrum dis- 
tortion of the CMB due to the Compton scattering from thermal electrons contained in the 
clusters of galaxies are obtained by solving the Kompaneets equation under the conditions 
that the probability of the multiple scattering of a photon is negligibly small and the electron 
temperature is much higher than the temperature of the CMB. This effect is known as the 
thermal Sunyaev-Zel'dovich effect (hereafter thermal SZE) (Sunyaev & Zel'dovich 1972). 

The spectrum distortion due to the plasma with a temperature gradient is deduced 
by the following way. In the plasma with a temperature gradient, the electron momentum 
distribution function is described as g(q) = g m (q) + Ag(q). Since the multiple scattering 
of a photon is negligible in the clusters of galaxies, the effects of the spectrum distortion of 
the CMB due to the Compton scattering from electrons belonging to g m (q) and electrons 
belonging to Ag(q) are decoupled. Therefore, each effects can be evaluated separately and 
the total spectrum distortion is given by the superposition of each results. Hereafter, the SZE 
due to electrons belonging to Ag(q) is referred to gradT SZE. Since the first non-zero moment 
of the Ag(q) appears from the third-order moment, to evaluate the effect of the Ag(q) 
on the CMB spectrum distortion the third order expansion in electron momentum of the 
collision term of the Boltzmann equation is required. The technical details are summarized 
in Appendix. A final result of third-order expansion is reduced to 



P) oi . °"t 

II • V ± e 
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(3) 



n . — - — ;n • VT e 

cut mc 6 



x 



X f (x, p) + f^r (x, p) + f^r (x, p) 



where x = ^ T , h is the Planck constant, c is a speed of light, v is the frequency of the 
CMB photon, T 7 is the temperature of the CMB, <tt is the Thomson cross section and n = £ 
is the direction cosine of the observed photon, respectively. Although the resultant photon 
distribution is anisotropic, the isotropic photon distribution for incident photons is assumed 
in the course of calculation since the multiple scattering is negligible and the incident photons 
are CMB photons. The quantities k B T e /mc 2 and E(p)/mc 2 are treated as the order of (v/c) 2 
variables. 



3. The spectrum of the gradient-T SZE 

The spectrum of the gradT SZE is obtained by inserting the Planck function in the right 
hand side of the above obtained third-order Boltzmann equation and integrating over line 
of sight. The frequency dependent amplitude of the distortion expressed in temperature is 
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obtained as 
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where d£ is the line of sight integral. Figure 1 compares the spectrum shape of the gradT 
SZE with the thermal and kinetic SZE(Hu et al. 1994; Birkinshaw 1999). Figure 2 shows 
the frequency dependence of the kernel function A VT (x) to see diagnostics of the spec- 
trum shape of the gradT SZE compared with the thermal and kinetic SZE(Hu et al. 1994; 
Birkinshaw 1999). The spectrum shows distinct characteristics compared with others. There- 
fore, the functional form of the non-Maxwellian part can be measured by this spectrum 
shape of the gradT SZE. It has a finite value at the cross over frequency of the thermal 
SZE, that is 218GHz. The location of maxima, cross over frequency, and minima are in- 
dependent of T e . They are at x max = 3.44 [y = 195GHz), x zero = 5.75; {y = 326GHz), 
and rr min = 8.36; [y = 437GHz), respectively. Since the amplitude of the spectral distortion 
is proportional to the inner product of the n with VT e , the sign of the spectrum depends 
on the relative direction of the line-of-sight to the temperature gradient. When the hotter 
region locates closer to the observer, the amplitude of the distorted spectrum is larger than 
the CMB in the Rayleigh- Jeans part. A spherically symmetric system does not show the 
spectrum distortion due to the gradT SZE since both of increment and decrement with the 
same amount appear in a single line-of-sight and these are cancelled out. When the line of 
sight is exactly perpendicular to the temperature gradient, there is no gradT SZE effect. 

The amplitude of the spectrum distortion due to the gradT SZE at the peak fre- 
quency is evaluated by the following way. The line of sight integral in the right hand 
side of the equation is reduced to It is independent from the electron density be- 

cause the heat conductivity does not depend on the electron density explicitly unless it is 
saturated (Cowie & McKee 1977). As decreasing the electron density the Coulomb mean free 
path becomes longer. In the region where the electron density is lower than 10" 5 cm -3 with 
/cbT 6 = lOkeV, the mean free path is longer than 2Mpc(Sarazin 1988). For the typical rich 
cluster, the electron density gets lower than 10~ 5 cm -3 in the region of 2Mpc apart from 
the cluster center. Therefore, in such a cluster out skirt the Coulomb mean free path is 
longer than temperature gradient scales and the heat conduction is saturated. As modeled 
by the Cowie and Mackee (1977), the saturated heat conductivity could be evaluated by the 
collisionless diffusion of the hot electrons with the Maxwellian velocity distribution. It is 
reasonable to assume that the electron velocity distribution is the Maxwellian in the cluster 
out skirt region where the heat conduction is saturated since the deviation of the velocity 
distribution function from the thermal equilibrium form is hard to be established without 
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any collision between particles. Therefore, we assume that Ag = in the cluster out skirt 
region where the heat conduction is saturated. For simplicity, we assume that Ag disappears 
suddenly at the points where the electron Coulomb mean free path becomes as long as the 
temperature gradient scale. The temperature of these transition points for the closer side 
and far side to the observer along the line-of-sight are denoted by T e 2 and T e i, respectively. 
Then, the amplitude of the gradT SZE is obtained as 



The result shows that the amplitude of the gradT SZE is proportional to the heat conduc- 
tivity. To extract the heat conductivity from the observed amplitude of the gradT SZE, the 
electron temperature in the cluster out skirt at the both side of the line-of-sight have to be 
measured. Since it is practically difficult, some assumptions must be made to extract the 
heat conductivity from the observed temperature map. The expected amplitude is two order 
of magnitude less than the thermal SZE and an order of magnitude less than the kinetic 
SZE but order of magnitude larger than the gravitational scattering effect. To distinguish 
the gradT SZE signal from other SZE signals, measurement of the spectrum by broad band 
or multi-band observations are required. 



The applications of the gradT SZE to three different situations in clusters of galaxies are 
considered in this section. First case is that the plasma has a smooth distribution of density 
and temperature. Second and third cases contain one discontinuity, such as a shock and 
a cold front along the line-of-sight, respectively. An idealized line-of-sight passing through 
clusters of galaxies is illustrated in Figure 3. The region denoted by D is a discontinuity. The 
regions denoted by 1 and 2 are smooth regions. The outer edges of these regions are defined 
by the place where the heat conduction starts to be saturated (Cowie & McKee 1977; Sarazin 
1988). For simplicity, the reduction rate of the heat conduction in these smooth regions, f K , 
is assumed to be constant through out the regions and suddenly becomes zero at each edges. 
The electron temperature at the edge of the discontinuity and at the outer edge of each 
smooth regions are denoted by suffixes D and o, respectively. Since the scale lengths of the 
discontinuities, L D , are comparable to or shorter than the Coulomb mean free path, the 
reduction rate of the heat conduction in the discontinuities, f' K , may be different from that 
in the smooth region and could be much less than 1. 





4. Applications to clusters of galaxies 



The general expression for the amplitude of the gradT SZE for the observer illustrated 



- 8- 



in Figure 3 is obtained by performing the line-of-sight integral as 



ATj r oc f K [T:W-T e D (2y/'} 

+f , [T D {2 y/2 _ t d (1)7/2] _ u [T o {1 y/2 _ t d (1)7/2] _ 



No discontinuity along a line of sight 

When there is no discontinuity along a line of sight, the reduction rate of the heat 
conductivity is constant through out the line of sight. By setting f' K = f K , the amplitude of 
the gradT SZE at v = 195GHz is reduced to 



The amplitude depends on the temperature at both outer edges of the smooth region and the 
heat conductivity in the smooth region. The amplitude is insensitive to any local variation 
of the temperature along line-of-sight. When the temperatures at the both outer edges are 
close to each other, the amplitude of the gradT SZE becomes zero. 

Across a shock front 

Although the electron distribution function in the cluster shock waves are quit uncertain, 
it is likely that the shock is collisionless and the thickness of the shock front is much smaller 
than the electron mean free path. Therefore, it is expected that the heat conductivity in the 
shock waves is much less than that in the smooth region, e.g. f' K <C f K . The amplitude of 
the gradT SZE across a shock front at v — 195GHz is reduced to 



where the effect comes from the difference of the electron temperatures at the both sides of 
the outer edges in the smooth regions is neglected. The adopted temperature of 30keV as 
the post shock gas temperature is somewhat larger than the average X-ray temperature of 
the cluster hot gas, that is 2-10keV. However, there are several observational and theoretical 
evidences which indicate that the temperature of the shock heated gas due to a cluster-cluster 
major merger could be as hot as 30keV. A cluster 1E0657-56 is an on-going merging cluster 
in which a contact discontinuity and a shock are identified (Markevitch et al. 2002). The 
temperature map of this cluster indicates the existence of the high temperature region hotter 
than 20keV. The combined analysis of the multi wave band SZE mapping observations of 
the most X-ray luminous cluster RXJ1347-1145 (Kitayama et al. 2004) showed that the 
inferred temperature of the south-east part of the cluster is well in excess of 20keV which is 
confirmed by the Chandra X-ray temperature measurement (Allen et al. 2002). A numerical 
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simulation of the head on collision of rich clusters (Takizawa 1999) showed that the shock 
heated gas temperature may temporary exceed 20keV. Therefore, in such an ideal case, the 
amplitude of the gradT SZE could be comparable to that of the kinetic SZE. An important 
point we have to emphasize is that the amplitude of the gradT SZE across the shock front is 
proportional to the heat conductivity in the smooth region. Therefore, this observation can 
be used to measure the heat conductivity in the smooth region. 

Across a cold front 

The thickness of cold fronts are yet unresolved. The current observed upper limits are 
comparable to the electron mean free path. Since the cold fronts are found in relatively 
large fraction of clusters, the cold fronts should not be a short lived transient structure. 
To maintain the cold fronts, the heat conductivity should be suppressed in many order of 
magnitude. It is likely that the heat conductivity in cold fronts is much less than unity, 
f' K <C 1, and much less than that in the smooth region, f' K <C f K . Therefore, the amplitude 
of the gradT SZE when the line of sight acrosses one cold front, is estimated as 



where the effect comes from the difference of the electron temperatures at the both sides of 
the outer edges in the smooth region is neglected. The dependence on the physical variables 
are the same as the case acrossing a shock. However, the temperature increase acrossing the 
cold front is not so significant since there is no thermalization of the kinetic energy as in 
the case of the shock. Therefore, the expected amplitude is much smaller than the case of 
acrossing the shock. 



We have shown that the inverse Compton scattering of the CMB photon with electrons 
in the ICM which has a temperature gradient, results in a new type of the spectrum dis- 
tortion of the CMB by the third-order perturbation theory of the Compton scattering. The 
spectrum has an universal shape. The cross over frequency appears at 326GHz. The sign 
of the spectrum depends on the relative direction of the line-of-sight to the direction of the 
temperature gradient. When the hotter region locates closer to the observer, the intensity 
becomes brighter than the CMB spectrum in the frequency regions lower than the cross over 
frequency. Therefore, the effect is distinguishable from other SZE effects by broad band or 
multi-bands observations of the SZE. The precise measurement of the spectrum shape of the 
gradT SZE provides an unique opportunity to measure the non-Maxwellian part of the elec- 
tron velocity distribution function when the ICM has a temperature gradient. The amplitude 




5. 



Discussion 
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of the spectrum distortion is proportional to the heat conductivity. It provides a possibility 
for a direct measurement of the heat conductivity. The expected amplitude of the gradT 
SZE signal is more than two order of magnitude smaller than the thermal SZE. The biggest 
signal for gradT SZE occurs when the line-of-sight across a shock front. The expected signal 
for this case is as large as ~ 30/iK when f K — 1 at maximum. This is ~ 10% of the thermal 
SZE. However, there are variety of the effects which distort the thermal SZE spectra, such 
as the relativistic correction (Birkinshaw 1999), the existence of the non-thermal electrons 
(Birkinshaw 1999), the hydrodynamical effect of the cluster merger (Koch 2004) and the 
turbulent motion giving rised by merger shock (Nagai, Kravtsov & Kosowsky 2003). These 
can become problematic when using multi-frequency separation techniques. 

The gradT SZE does not cause any polarization since the quadrapole moment of the 
CMB intensity distribution averaged over all the electron momentum is zero. The linear po- 
larization of the Compton scattered photons is generated by the systematic transverse motion 
of the electron system since it leads a finite quadrapole moment of the CMB intensity distri- 
bution in the frame moving with the same velocity as the transverse motion of the electron 
system(Sunyaev & Zel'dovich 1980). As explained in the Section 1 and appendix, the first 
and second moments of the electron momentum on Ag(q) are zero(Okabe & Hattori 2003). 
Therefore, the Ag(q) does not generate a quadrapole moment of the CMB intensity distri- 
bution. 

The third order Doppler effect could provide the physical explanation of the spectrum 
distortion in the gradT SZE although it is not satisfactory enough. The scattering can be 
treated as almost elastic in the electron rest frame, i.e., p' = p. The transformation into the 
laboratory frame induces a Doppler shift and an energy transfer to the photon amount of 5p 
given as 

5p 1-P-n 

P 1-P'U' 

Averaging over the incoming photon direction n which is isotropically distributed, we obtain 

The averaging over electron velocity with velocity distribution function of Ag results in 

Sp\ 6 n-VT e f\ e \v th Ak B T 
f K K S p ~ ±0.14 



p I 5n e m e c 3 \ L J c m e c 2 ' 

where A e is the Coulomb mean free path of electron, L is the temperature gradient scale 
length and v th = ^/2k B T/m e is the electron thermal velocity. This is an averaged fractional 
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energy shift of the photon by a single scattering. When the hot region locates foreground 
of the cold region along the line of sight, the minus sign in the last equality is taken. The 
scattered photons lose their energy. Therefore, the original Planck spectrum shifts toward 
low frequency region contrary to the case of the thermal SZE. This results in the increment 
in the Rayleigh- Jeans part and the decrement in the Wien part contrary to the thermal SZE. 
When the cold region locates foreground of the hot region, the plus sign in the last equality 
is taken. The scattered photons gain the energy. In this case, the original Planck spectrum 
shifts toward high frequency region. The resultant spectrum distortion is the decrement in 
the Rayleight- Jeans part and the increment in the Wien part as same as the case of thermal 
SZE. 
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Appendix 

Some technical details for performing the third order perturbation theory for the Comp- 
ton scattering are summarized. In the following discussion, p = E(p) and p' = E(p') are 
used for simplicity. 

The third order expansions of each terms contained in the collision integral are obtained 
as follows. The expansion to third order of the energy conservation can be treated as an 
expansion in Ap = E(q) — E(q') as 

S(p + q—p' — q') = 5(p'-p-Ap) 

d 

= S(p-p)- Ap—5(p - p) 

Id 2 1 d 3 

+ 2 Ap2 dp^ 6 ( p> ~ P) ~ 6 Ap3 dp^^ P ' ~ P) ' 
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Since the electron momentum distribution function does not depends on q in the first term 
of the collision integral, first performing the integral by q and q is replaced with q' + p' — p 
by using the momentum conservation. Then, Ap is written as 

Ap = swyT7 (p ' (p '- p)+q, ' (p '- p)) - 

For the second term in the collision integral, the integral by q' is performed first. The energy 
difference Ap is written as 

In the pre-scattering electron rest frame, the energy of the scattered photon is described by 
the energy of the incident electron as 

P 

P l + ^l-n-n')- 

The Doppler formula provides a relation between the photon energy in the rest frame and 
that in the laboratory frame as 

p = 7 p(l-/3-n), p' = 1P '(l-f3 -n'), 

where 7 is the Lorentz factor and (3 = q/mc is the incident electron velocity in the laboratory 
frame relative to the speed of light. The Lorentz invariant relation of p-p' = p-pf provides a 
relation of the photon scattering angle between the rest frame and the laboratory frame as 

_ ~ -/3-(n + n') + (/3-n)(/3-n')+/3 2 (l-n-n') + n-n' 

l-/3-(n + n') + (/3-n)(/3-n') ' 

Since the ratio of the photon energy to the electron rest mass energy is order of [3 2 , the third 
order expansion of \M\ in (5 is expressed as 

\M\ 2 = 2(4vr) 2 e 2 (l + fi-n') 

= 2(4tt)V[1 + (n • n') 2 - 2/3 • (n + n')n ■ n'(l - n • n') 

+(/3 • (n + n')) 2 (l - 3n • n')(l - n • n') + 2/3 2 n • n'(l - n ■ n') 

+2(f3 -n)(f3- n')n ■ n'(l - n ■ n) + 2(/3 ■ (n + n)) 3 (l - 2n • n')(l - n ■ n) 

-2((3 ■ (n + n')),3 2 (l - n • n')(l - 2n ■ n) 

-2((3 ■ (n + n'))(/3 • n)(/3 • n')(l - 3n • n')(l - n ■ n')]. 



The third order expansion of the electron energy is E(q) ~ mc 2 (l + q 2 /2m 2 c 2 ). 
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The 1st and 2nd moments of electron momentum for the Ag are zero, and the 3rd 
moments are obtained as 

--m 2 f K K S pd x T e , for % = x, j = k = y or z, 
o 

6 

--m 2 f K K SP d x T e , for i= j = k = x 
5 

0, otherwise 

where the x axis is taken to the direction of the temperature gradient. The no n- zero terms 
in the third order expansion of the collision term are only the third order moment of electron 
momentum and can be evaluated by the above results. The evaluations of the each terms 
are straight forward. 



/ 



d 3 q 

J2nf 
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Fig. 1. — Comparison of the CMB spectrum distortion due to the gradT SZE with the 
thermal SZE and the kinetic SZE. The absolute amplitude of each effects are compared 
by log-linear plot. A vertical axis is described in K. A horizontal axis is the frequency of 
the CMB photons normalized by its temperature. The solid line shows the spectrum of 
the gradT SZE with amplitude shown in Eq.(6). The dashed line shows the spectrum of the 
thermal SZE assuming an uniform sphere of radius 0.5Mpc with electron density of 10 _3 cm -3 
and electron temperature of lOkeV. The dotted line shows the spectrum of the kinetic SZE 
assuming the same uniform sphere moving with a peculiar velocity of 600km/s. 
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x=hv/k B T y 

Fig. 2. — The kernel functions of the three SZE spectrums are shown. The solid line 
shows the kernel function of the gradT SZE. Dashed line and dotted line are the kernel 
function of the thermal SZE, A T (x) = x A e x /(e x — l) 2 (xcoth(a;/2) — 4), and the kinetic 
SZE, A (x) = x 4 e x /(e x — l) 2 , respectively(Birkinshaw 1999). A vertical axis is described in 
dimensionless arbitrary unit. 
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Fig. 3. — An illustration of the line of sight which passes one discontinuity. The line of 
sight toward an observer is illustrated by an long arrow. The region denoted by D is a 
discontinuity. The regions denoted by 1 and 2 are the regions where the temperature and 
density distributions of the hot gas are smooth. The electron temperatures at the edges 
of the discontinuity is denoted by the super script D and the electron temperature at the 
outer edges of the cluster is denoted by the super script o, respectively. The fractional heat 
conductivity in the smooth region is described by f K and that in the discontinuity is described 
by f' K , respectively. 



